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Abstract: We compute the rate for diffractive T meson production at the Tevatron and 
the LHC. The T is produced diffractively via the subprocess j+p — > T+p where the initial 
photon is radiated off an incoming proton (or antiproton) . We consider the possibility to 
use low angle proton detectors to make a measurement of the jp cross-section and conclude 
that a measurement of the cross-section at a centre of mass energy in excess of 1 TeV is 
possible at the LHC. This is in the region where saturation effects are likely to reveal 
themselves. 
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1. Introduction 

In this paper we calculate the rate for the reaction p+p — > p + T{nS) +p where n = 1, 2, 3. 
The T is produced diffractively after one of the incoming protons radiates a photon, as 
illustrated in Figure 1, i.e. via the sub-process 7 +p — > T + p. 

The reaction is interesting primarily because the LHC will probe values of the 7]? 
centre-of-mass energy well beyond those that were reached at HERA. As a result, there 
is an opportunity to examine QCD in the region where non-linear (saturation) effects are 
expected to be important. In addition, the odderon is expected to contribute at some level 
and an accurate measurement would help to constrain the theory. We shall not consider 
the odderon in this paper. 

Ordinarily this is a measurement that could only be performed at low LHC luminosities 
(~ 10 31 cm _2 s _1 ) in order to avoid problems with pileup. However, there is interest 
in supplementing the LHC general purpose detectors with low angle proton detectors at 
420 m from the interaction points; the primary goal being the study of 'central exclusive' 
production of new physics, and in particular the reaction p+p — > p+H+p [1]. The T states 
are too light to be produced in conjunction with two measured protons but measurement 
of one proton is feasible and that may be sufficient to control pileup induced backgrounds 
at instantaneous luminosities in excess of (10 33 cm _2 s _1 ). Work is currently in progress 
to measure exclusive T production at the Tevatron [2] and the CDF collaboration has 
identified centrally produced charm meson states [3] and centrally produced dileptons [4] . 

There already exist a number of papers predicting the rate for this process [5, 6, 7, 
8, 9, 10]. We add to what has been done by providing predictions using the dipole model 
cross-section of [11] that successfully describes a wide range of the HERA data [12]. We 
also make an estimate of the effect of detector acceptances and the impact of measuring 
one of the protons. 

In the following section we explain how to calculate the rate for this process. The 
flux of photons off an incoming proton is well known, as is the dipole cross-section that 
determines how the photon produces the T after first fluctuating into a qq pair. The latter 
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Figure 1: Diffractive photoproduction of the T meson in pp collisions. At the Tevatron one p 
stands for an antiproton and the other for a proton while at the LHC, both stand for protons. The 
T decays into a fi + p~ pair which is detected. At the LHC, there is opportunity to also detect one 
of the protons. 



is well determined , primarily as a result of high quality data on the deep inelastic structure 
function F2 measured at HERA [13, 14]. In Section 3, we compute the photoproduction 
cross-section and compare our predictions with the available HERA data. In Section 4, we 
present our predictions for the rapidity distributions of the T and the expected rates at 
the Tevatron and the LHC, including the effects of proton tagging. 



2. The pp cross-section 

Entirely in analogy to photoproduction at HERA, where almost on-shell photons are radi- 
ated off incoming electrons according to the Weisz acker- Williams distribution, protons at 
the LHC can radiate almost real photons that can then scatter off protons heading in the 
opposite direction. A significant fraction of these jp reactions will be diffractive and leave 
the incoming proton intact. Amongst these will be the process 7 + p — > V + p where V is 
a vector meson and that is where our interest lies. 

If the T is produced at rapidity Y then its energy Ey « My cosh Y (we neglect the 
transverse momentum of the meson) . Defining 

<-^. y (2.1) 
where y/s is the center-of-mass energy of the pp system, the cross-section of interest is 

dM dY^ TP) = tfr/^' Q2) + (Y^-Y), (2.2) 

where / 7 / p (£, Q 2 ) is the photon flux given by [15]: 

f U ^_"cm l + (l-Q 2 1 1 , . 



x At least for light quark scattering. 
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and /i 2 = 0.71 GeV 2 fixes the electromagnetic form factor of the proton. £ is the fractional 
energy loss of the proton. The cross-section for 7+p^T+pis denoted by a^ p {W 2 ) where 
W is the 7p centre of mass energy, i.e W 2 = M-f\/sexp(±y). W is not uniquely defined 
by the T's rapidity since any one of the two incoming protons can radiate the photon. The 
two terms on the right-hand-side of Eq. (2.2) allow for both possibilities. 2 
The minimum virtuality of the photon is 

Qmin ~ t 2 ™ 2 p , (2.4) 
where m p is the proton mass. Integrating Eq. (2.3) over Q 2 yields the integrated flux [15]: 
« cm 1 + (1 - Q 2 ( 11 , 3 3 1 \ 

with 

A(0 = l + V 2 /QLn (2-6) 
and thus the rapidity distribution of the T is given by 
da(pp — ► pTp) 



dY 



e/ 7 / P (0 v-rpfa) + (Y^-Y). (2.7) 



So much for the photon flux, we are now ready to tackle the cross-section for 7 + 
p — > T + p. If the invariant mass W 2 3> My then we are in the diffractive regime and 
consequently the imaginary part of the forward (t = 0) amplitude can be approximated 
by 3 : 

9m A 1P (W 2 ) = ±£ / dV d2 ^xk^ r K,Hh^ r M*, r ) > ( 2 - 8 ) 
X,h,h 

where ^ ^("^ r ) an d r ) are ^ e light-cone wavefunctions of the photon and the T 

(of helicity A). They represent the amplitude for the vector meson to fluctuate into a bb 
pair of transverse size r and with the quark carrying energy fraction z and helicity h = ± ^ . 
We take x = M 2 /W 2 . 

The cross-section cr(x, r) is called the dipole cross-section and it represents the cross- 
section for the bb pair to scatter off a proton. The formalism we are describing has been 
very successful in explaining a very wide range of HERA data, including the extremely 
precise data on the Fi structure function at low x [12]. The dipole cross-section has 
been extracted, for light-quark dipoles, from HERA data and we use a parameterization 
introduced in [11]; we shall turn to this shortly. First we address the calculation of the 
light- cone wavefunctions. 

The photon wavefunction is well known and, since the longitudinal wavefunction of 
a real photon vanishes, we shall only need the transverse (A = ±1) photon wavefunction 
given by 

- i(2h)S h _- h e iX * [(1 - z)5 x ^ 2h + z5 x ,2h] m^m^)} , (2.9) 



2 By adding cross-sections, we are neglecting the interference between the two amplitudes. 
3 We assume that the T retains the helicity of the photon. 
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where e 2 = 47ra em and ej = —1/3 is the electric charge of the b quark. 

The vector meson wavefunction is less well known. Various models are discussed in 
[16] and here we shall use the boosted gaussian wavefunction, which works well for the 
light mesons, p and <j> [16], and also for the heavier J/ty meson [12]. In this model, the 
meson light-cone wavefunction is given by 



+ i(2h)S h> _- h e^ [(1 - z)S x ^ 2h + z5 x>2h ] d r }4> n s(z, r) 



(2.10) 



where (f) n s(z, r ) is obtained by boosting a Schrodinger gaussian wavefunction using the 
Brodsky-Huang-Lepage prescription [17]: 



0sch(p 2 ) -»■ 4>nS 



k 2 + m 2 

4z(l - z) 



m, 



(2.11) 



and p is the relative three-momentum between the quark and antiquark in the meson's rest 
frame while k is their relative transverse momentum in the boosted frame. The subscript 
nS reminds us we are to treat the different T states seperately. After a two-dimensional 
Fourier transformation of the resulting light-cone wavefunction to transverse r-space, we 
obtain 

"n-l 

J2^nS,kRls D 2k (r,z) 
Jc=0 



4>n,s{r, z) 
with a n so = 1- The operator 



G n s(r, z) 



m 2 f -Vl 
D(r, z) = -y^- - mj 



4z(l - z 

with V 2 = ^d r + d 2 , acts on the Gaussian function 

mlRl c \ ( 2z (l - z)r 2 
R ls 



G nS (r, z) = M nS z(l - z) exp ( 



exp - 



exp 



Explicitly we obtain 4 



4>is{r,z) = Gi S (r,z) , 



(2.12) 



(2.13) 



(2.14) 



(2.15) 



<f>2s(r, z) = G 2 s(r, z)[l + a 2 s,i 92s(r, z)] 



and 



4>3s(r, z) = G 3S (r, z) \ 1 + 035,1 g 3S (r, z) + a 3S ,2 
where 



9is(r,z) + 4:[l 



4z(l - z)r' 

^35 



2r)2 t m}R 2 nS Az{l- Z y 



g nS (r,z)=2-m f R nS + —_ z) r2 



nS 



(2.16) 



(2.17) 
(2.18) 
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Figure 2: The scalar part of the light-cone wavefunction for T(15) (left), T(25) (center) and 
T(35) (right) as a function of the transverse size, r (GeV -1 ), of the bb pair with the light-cone 
momentum fraction carried by the quark, z = 0.5 (solid), z = 0.3 (dashed) and z = 0.2 (dotted). 

The scalar wavefunctions <j)\Si 4>is an d 4>3S are shown in Figure 2 for z = 0.2, 0.3 and 0.5. 

The parameters a n s t k (f° r k > 0) have been fixed by requiring that the wavefunctions of 
different states be orthogonal to each other. The values of R n s are fixed so as to reproduce 
the experimentally measured electronic decay widths [16]. The values of the parameters 
thus obtained, together with the predicted decay widths, are given in Table 1 and the 
resulting wavefunctions are shown in Figure 3. 

Boosted Gaussian parameters 



n 


R ls 




"nS,2 




r e + e - 


■pexp 

e+e - 


1 


0.567 






0.481 


1.340 


1.340 ±0.018 


2 


0.831 


-0.555 




0.624 


0.611 


0.612 ±0.011 


3 


1.028 


-1.219 


0.217 


0.668 


0.443 


0.443 ± 0.008 



Table 1: Parameters and predicted electronic decay widths of the Boosted Gaussian light-cone 
wavefunctions in appropriate GeV based units. The decay widths are given in kcV. Experimental 
values are taken from [19] and we take mi, = 4.2 GeV. 

It remains to specify the parametric form of the dipole cross-section. Here we use the 
saturation model presented in [11] (which we subsequently refer to as the "FSSat" model). 
It has the following form: 

a{x, r) = Aht 2 x~ Xh for r < ro and 

= A s x~ Xs for r > n. (2.19) 

In the intermediate region ro < r < r\, the dipole cross-section is determined by interpo- 
lating linearly between the two forms of Eq. (2.19). ro is fixed to be the value at which the 
hard component is some specified fraction of the soft component, i.e. 

(Hp, r )/a{x, n) = f, (2.20) 

4 We note that if the double partial derivative <9j? is used in Eq. (2.13) instead of the Laplacian operator, 
V^, the 25 wavefunction reduces to that presented in [18]. 
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Figure 3: The light-cone wavefunction squared, I 1 ? 1 "! 2 , for the transversely polarised T(15) (left), 
T(25) (center) and T(35) (right) as a function of the transverse size, r (GeV -1 ), of the bb pair and 
the light-cone momentum fraction, z, carried by the quark. 



where / is treated as a parameter that is determined, along with the other parameters 
As, Ah, Ag, Xh and n, by a fit to the total structure function F2 data from HERA [11]. 

We now have almost all of the ingredients we need in order to compute the amplitude 
(2.8) for T production. The forward differential cross-section is given by 



da. 



dt 



t=o 



1 

167T 1 



\A 



(2.21) 



and the total cross-section is obtained by assuming an exponential fall-off with increasing 
\t\, i.e 

1 dd^ri 

(2.22) 



a. 



B r dt 

where the diffractive slope, By, is taken to be [16] 

14.0 



t=o 



B r = N 



(M T /GeV) - 4 



+ 11= 3.68 GeV" 2 



(2.23) 



with N = 0.55 GeV~ 2 . Before presenting our predictions for hadron-hadron collisions, we 
shall first compute the photon-proton cross-section a lp and compare it to existing HERA 
data. 



3. The 7p cross-section 

To compare to the data, we compute 

3 

°7P = a ip->~£(nS)p x B n , (3.1) 
n=l 

where B n is the branching ratio of the T(nS) to muons. In Figure 4, we compare our 
predictions to the HERA data [20, 21, 22]. 

As can be seen, the theory curve (dotted) lies significantly below the HERA data. 
However, we have ignored two important corrections. The first is the contribution from 
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Figure 4: Predictions for the 7p cross-section. Dotted curve - no skewedness correction and 
without the real part; dashed curve - including the real part; dot-dashed curve - including the real 
part and a skewedness correction. The solid curve is the dot-dashed curve normalised to the HERA 
data and the double-dot-dashed curve is the parameterization to the HERA data described in the 
text [6]. 



the real part of the amplitude and including it increases the cross-section for each state by 
a factor of (1 + f3 2 ) where (3 is the ratio of the real to imaginary part of the jp amplitude, 
given by 

P = tan (|A) (3.2) 

with 

_ dln|5m A\ 
A " 31n(l/s) [6 ' 6) 

where A is given by Eqn. (2.8). 

Figure 5 shows the ratio of the real to imaginary parts of the amplitude for each state. 
We note that the calculation of this ratio is model-dependent especially for the higher 
states. The dashed curve in Figure 4 shows the effect of including the correction due to 
the real part. 

The second correction arises because the amplitude is not diagonal - the mass of the T 
is large and timelike compared to the small, spacelike, photon virtuality. In contrast, our 
dipole scattering cross-section was extracted from data at zero momentum transfer. We 
estimate the corrections from this source by multiplying the amplitude by a factor of [23] 

2 2A+3 F(A + 5/2) 

R ° {x) - ^ r(A + 4) • (3 - 4) 



Our prediction, with this correction included, is shown in Figure 4 as the dot-dashed curve. 

A few comments are in order regarding our estimation of the above two corrections. 
Both of them depend on A, which we estimate using Eqn. (3.3). We have checked that 
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Figure 5: The ratio of the real to imaginary part of the amplitude for the different states. 

our results do not vary very much upon choosing a fixed value of A = 0.3. Since the T 
wavefunction extends to larger r for the excited states, relative to that of the Is state (see 
Figure 3), one might naively expect that the energy dependence is softer for the former, 
resulting in a lower effective value of A. This is actually not the case as one can infer 
from Figure 5. This is because the excited state wavefunctions have nodes, which leads 
to a partial cancellation between the soft and hard parts of the amplitude which in turn 
results in an effective harder energy dependence (higher A) in the amplitudes. The degree 
of cancellation is rather sensitive to the model used for the dipole cross-section. This has 
been referred to previously as the node effect [18]. 

Even with these corrections, the theory curve still lies below the data. This is not a 
particularly surprising result and it is typical of other dipole model calculations [6, 24]. 
Quite possibly, one could do better by improving the vector meson wavefunctions and in 
that case the uncertainty would mainly influence the normalisation of the cross-section. 
Moreover, the value chosen for the b quark mass also affects the overall normalisation. For 
these reasons, we follow the authors of [6] and rescale our results (by a factor 2.0 with 
rrib = 4.2 GeV) in order to achieve optimum agreement with the data. The result is the 
solid curve in Figure 4. By doing this, we account also for the uncertainties in the real part 
and skewedness corrections in a purely phenomenological way. 

In [6], the parametrisation 



was found to fit the HERA data for Y(IS'). Both HI and ZEUS measure a* p and then 
extract the IS cross-section by assuming that the ratios of cross- section times branching 
ratio are the same as those measured by the CDF collaboration [25], i.e. (<72S • -B2s)/(cis ' 
Bis) = 0.281 ± 0.0484 and (a 3S ■ B 3S )/(a ls ■ B 1S ) = 0.155 ± 0.0319. We use Eq. (3.5), 
together with the CDF ratios, to produce the double-dot-dashed curve in Figure 4. 

We can also compute the cross-section ratios between the different T states in order to 
compare to the CDF data. These are shown in Figure 6. Our results are below the values 
obtained by CDF quoted above. For the 2S : IS" ratio, in the range 50 < W < 200 GeV, our 
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Figure 6: The ratio of cross-sections for the different states. 



prediction is also below the value calculated in [7] using perturbative QCD and a Gaussian 
wavefunction. However, theoretical predictions for these ratios are rather uncertain since 
the 2S : 15 ratio is very sensitive to the details of the meson wavefunction and also to the 
value of A. 

Note that the FSSat model predicts an energy dependence that is less steep than the 
W 1 '® dependence of the fit. We might take the difference to be indicative of the challenge 
facing future experiments, i.e. they should be able to distinguish between the two. 



4. Results 



We are now ready to present our predictions for the hadroproduction cross-section. We 
assume that the angular distribution of the decay muons in the T rest frame is oc (1 + 
cos 2 9* ) , which corresponds to a distribution of 



d 2 iV 3 l-(3 2 

d(cos(9)d</> ~ i67r(l - (3cos9) 2 



1 + 



cos 9 — (3 
(1-/3 cosd 



(4.1) 



in the lab: (3 = tanh Y is the speed of the T and 9 is the polar angle relative to the direction 
in which the T travels when (3 > 0. 5 Thus 



d 2 a(pp — > }p) dN 



d(cos#)dY 



d(cos 9) 



(4.2) 



We integrate over 9 such that both the n and /i + are emitted at an angle greater than 
#min relative to the beam. To do that we need a, the emission angle of the 



cos a 



(1 + f3 2 ) cos 9 -2(3 
1 + (3 2 -2(3 cos 9 ' 



(4.3) 



Recall that we neglect the transverse momentum of the meson, which means it travels along the beam 
axis. 
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Figure 7: The T rapidity distribution without any cuts on the final state particles at the Tevatron 
(left) and LHC (right) using the FSSat dipole model and the parameterization the photoproduction 
data described in the text. 

The transverse momentum of the muon (or anti-muon) is given by 



and we also cut on a minimum value of P™ m for both leptons. 

For the Tevatron, we approximate the angular acceptance of the CDF muon detector 
using 9 m [ n = 33.5°, whilst for the D0 detector we take 8 m i n = 15.4°. At the LHC we 
assume, 9 m \ n = 7.7° (see also [26]). In all cases we show results for p™ m = 3 GeV and 
P t mm = 4 GeV. At the LHC, we may also have the possibility to detect one of the protons 
in the region 420 m from the interaction point (there will not be any acceptance to detect 
both). According to the studies presented in [1], forward detectors in the 420 m region have 
acceptance for protons with fractional energy loss 0.002 < £ < 0.018 on the right-hand side 
and 0.0015 < £ < 0.014 on the left-hand side. 

Figure 7 shows the T rapidity distributions at the LHC and the Tevatron respectively 
before any cuts have been applied. There is a striking difference (especially at the LHC) 
between the shapes of the distributions obtained using the the FSSat dipole model and the 
parameterization of the photoproduction cross-section that is based on the HERA data. 
One could therefore hope that measurements of the rapidity distributions would be able 
to discriminate between dipole models and thereby contrain the gluon distribution. In 
Figure 8 we indicate the effect of cutting on the muon pt and rapidity and, in the case of 
the LHC, the effect of observing one of the protons in the proposed 420 m detectors. Our 
results show that although the strong sensitivity to the energy dependence on the dipole 
cross-section is diminished, it is still large enough that one could hope to constrain the 
theory. 

Tagging one of the protons does severely limit the acceptance of a measurement: Cen- 
trally produced mesons are cut out because the 420 m detectors select events in which the 
measured proton loses much more energy than the unmeasured proton. The boost is so 
large in fact that the upsilon always travels in the direction of the tagged proton. However, 
detecting a proton does have the potential advantage that the measurement need not only 



Pt 



Ml sin 9 



(4.4) 



2E T (1 - p cos 9) 
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Figure 8: The T rapidity distribution after cuts at the Tevatron (top) and LHC (bottom) using 
the FSSat dipole model and the the parameterization the photoproduction data described in the 
text. 



Cross-sections at the Tevatron 





0"FSSat(fb>) 


0>it(fb) 


CDF 


350 


406 


D0 


685 


838 



Table 2: The predicted cross-sections for the Tevatron using the FSSat dipole model (first column) 
and the power-law fit to the photoproduction data (second column). Cuts on the muon rapidities 
and transverse momentum have been applied (P t > 4 GeV). 



be performed using data collected at low (< 10 32 cm~ 2 s _1 ) luminosities. With one tagged 
proton, one could hope to control the pileup background after utilising the fact that the 
two charged muons should point to a single vertex (and no other tracks point to the same 
vertex) and that they should also combine to produce an T rapidity that is in agreement 
with the value inferred from the extracted (one measured directly and one inferred) proton 
momentum fractions. Further detailed simulations would be required to establish whether 
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Cross-sections at the LHC 



No. of tagged protons 


0"FSSat(ft>) 


ffPit(fb) 





5133 


10351 


1 


3035 


6802 



Table 3: The predicted cross-sections for the LHC using the FSSat dipole model (first column) 
and the power-law fit to the photoproduction data (second column). Cuts on the muon rapidities 
and transverse momentum (P t > 4 GeV) have been applied (first row). The result of additionally 
measuring one of the outgoing protons is shown in the second row. 

this method will work at 2 x 10 33 cm _2 s _1 and above. In any case, it is worth noting that 
10 32 cm _2 s _1 corresponds to 1 fb _1 per year. That translates to a total of over 5000 signal 
events in an environment where there will be, on average, less than one pp interaction per 
bunch crossing. This should be sufficient to produce a measurement that will constrain 
QCD models of saturation. 

Finally, it is worth pointing out that should it be possible to make a cut on the measured 
proton's transverse momentum then it would become possible to make a direct measure- 
ment of the photoproduction cross section. This is so since the transverse momentum of 
protons that radiate photons is typically much smaller than the transverse momentum of 
protons that do not, i.e. Eq. (2.3) is much softer than the exp(Bj;t) dependence implied by 
Eq. (2.23). Requiring that the transverse momentum of the tagged proton p? < 100 MeV 
is at least 60% efficient for selecting events in which the tagged proton radiates a photon 6 . 
The contamination from events in which the non-tagged proton radiates a photon would 
be only 4%. If the cut is lifted to 300 MeV, the ratio decreases to 89% : 28%, which still 
constitutes a very significant enhancement. 

After making such a cut, since we now have an enriched sample of events in which 
the tagged proton radiated the photon, we can eliminate the photon flux and extract the 
photoproduction cross-section. Figure 9 illustrates the possibilities. The region at larger 
W (much larger than can be probed at HERA and the Tevatron) is obtained by insisting 
that the tagged proton pr be sufficiently small, in which case the tagged proton most likely 
radiated the photon. The lower W region can be measured by making the complimentary 
cut, i.e. by insisting that the proton have pr above some value. This is equivalent to 
assuming that the untagged proton emitted the photon. In this way, the proton detectors 
facilitate a measurement of the jp — ► Tp cross-section at W = 1 TeV, which is well in the 
range where saturation effects are expected to reveal themselves. 

So far, we have made no mention of the issue of "gap survival" and as far as the 
overall rate is concerned it is not expected to provide much suppression, since the collision 
is typically rather peripheral. However, the gap survival should depend rather strongly 
on the proton pr, with gaps being filled in more often at larger values of pr- Indeed 
measuring the px distribution has been advertised as a means to probe the gap-filling 
mechanism. This physics would need to be under control before one could reliably extract 

6 The efficiency depends weakly upon the energy lost by the proton through Eq. (2.4). 
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Figure 9: The photoproduction cross-section extracted after cutting on the transverse momentum 
of the tagged proton. 

the photoproduction cross-section by exploiting a pt cut on the scattered proton [27, 28] . 
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